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Abstract 

An extension of the Standard Model by three singlet fermions with masses smaller than the 
electroweak scale allows to explain simultaneously neutrino oscillations, dark matter and baryon 
asymmetry of the Universe. We discuss the properties of neutral leptons in this model and the 
ways they can be searched for in particle physics experiments. We establish, in particular, a lower 
and an upper bound on the strength of interaction of neutral leptons coming from cosmological 
considerations and from the data on neutrino oscillations. We analyse the production of neutral 
leptons in the decays of different mesons and in pp collisions. We study in detail decays of neutral 
leptons and establish a lower bound on their mass coming from existing experimental data and Big 
Bang Nucleosynthesis. We argue that the search for a specific missing energy signal in kaon decays 
would allow to strengthen considerably the bounds on neutral fermion couplings and to find or 
definitely exclude them below the kaon threshold. To enter into cosmologically interesting param- 
eter range for masses above kaon mass the dedicated searches similar to CERN PS191 experiment 
would be needed with the use of intensive proton beams. We argue that the use of CNGS, NuMI, 
T2K or NuTeV beams could allow to search for singlet leptons below charm in a large portion of 
the parameter space of the z^MSM. The search of singlet fermions in the mass interval 2 — 5 GeV 
would require a considerable increase of the intensity of proton accelerators or the detailed analysis 
of kinematics of more than 10^° B -meson decays. 



PACS numbers: 14.60.Pq, 98.80.Cq, 95.35.-hd 



I. INTRODUCTION 



In a search for physics beyond the Standard Model (SM) one can use different types 
of guidehnes. A possible strategy is to attempt to explain the phenomena that cannot be 
fit to the SM by minimal means, that is by introducing the smallest possible number of 
new particles without adding any new physical principles (such as supersymmetry or extra 
dimensions) or new energy scales (like the Grand Unified scale). An example of such a theory 
is the renormalizable extension of the SM, the z/MSM (neutrino Minimal Standard Model) 
[1, 2], where three light singlet right-handed fermions (we will be using also the names 
neutral fermions, or heavy leptons, or sterile neutrinos interchangeably) are introduced. 
The leptonic sector of the theory has the same structure as the quark sector, i.e. every left- 
handed fermion has its right-handed counterpart. This model is consistent with the data 
on neutrino oscillations, provides a candidate for dark matter particle - the lightest singlet 
fermion (sterile neutrino), and can explain the baryon asymmetry of the Universe [2]. A 
further extension of this model by a light singlet scalar field allows to have inflation in the 
Early Universe [3]. 

A crucial feature of this theory is the relatively small mass scale of the new neutral 
leptonic states, which opens a possibility for a direct search of these particles. Let us review 
shortly the physical applications of the z/MSM. 

1. Neutrino masses and oscillations. The z/MSM contains 18 new parameters in compar- 
ison with SM. They are: 3 Majorana masses for singlet fermions, 3 Dirac masses associated 
with the mixing between left-handed and right-handed neutrinos, 6 mixing angles and 6 
CP- violating phases. These parameters can describe any pattern (and in particular the 
observed one) of masses and mixings of active neutrinos, which is characterized by 9 param- 
eters only (3 active neutrino masses, 3 mixing angles, and 3 CP- violating phases). Inspite of 
this freedom, the absolute scale of active neutrino masses can be established in the z/MSM 
from cosmology and astrophysics of dark matter particles [1, 4-7]: one of the active neu- 
trinos must have a mass smaller than O{10~^) eV. The choice of the small mass scale for 
singlet fermions leads to the small values of the Yukawa couphng constants, at the level 
10~^ — 10~^^, which is crucial for explanation of dark matter and baryon asymmetry of the 
Universe. 

2. Dark matter. Though the z/MSM does not have any extra stable particle in comparison 
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with the SM, the hghtest singlet fermion, A^i, may have a hfe-time tn^ greatly exceeding the 
age of the Universe and thus play a role of a dark matter particle [8-11]. Dark matter sterile 
neutrino is likely to have a mass in the 0(10) keV region. The arguments leading to the 
keV mass for dark matter neutrino are related to structure formation and to the problems 
of missing satellites and cuspy profiles in the Cold Dark Matter cosmological models [12- 
15]; the keV scale is also favoured by the cosmological considerations of the production of 
dark matter due to transitions between active and sterile neutrinos [8, 9]; warm DM may 
help to solve the problem of galactic angular momentum [16]. However, no upper limit on 
the mass of sterile neutrino exists [3, 17] as this particle can be produced in interactions 
beyond the z/MSM. The radiative decays of A^^i can be potentially observed in different X-ray 
observations [10, 18], and the stringent limits on the strength of their interaction with active 
neutrinos [5, 19-27] and their free streaming length at the onset of cosmological structure 
formation [28-31] already exist. An astrophysical lower bound on their mass is 0.3 keV, 
following from the analysis of the rotational curves of dwarf spheroidal galaxies [32-34]. 
The dark matter sterile neutrino can be searched for in particle physics experiments by 
detailed analysis of the kinematics of /5 decays of different isotopes [35] and may also have 
interesting astrophysical applications [36]. 

3. Baryon asymmetry of the Universe. The baryon (B) and lepton (L) numbers are not 
conserved in the i/MSM. The lepton number is violated by the Majorana neutrino masses, 
while B + L is broken by the electroweak anomaly. As a result, the sphaleron processes 
with baryon number non-conservation [37] are in thermal equilibrium for temperatures 100 
GeV <T < 10^^ GeV. As for CP-breaking, the i/MSM contains 6 CP-violating phases in 
the lepton sector and a Kobayashi-Maskawa phase in the quark sector. This makes two 
of the Sakharov conditions [38] for baryogenesis satisfied. Similarly to the SM, this theory 
does not have an electroweak phase transition with allowed values for the Higgs mass [39], 
making impossible the electroweak baryogenesis, associated with the non-equilibrium bubble 
expansion. However, the i/MSM contains extra degrees of freedom - sterile neutrinos - which 
may be out of thermal equilibrium exactly because their Yukawa couplings to ordinary 
fermions are very small. The latter fact is a key point for the baryogenesis in the i/MSM, 
ensuring the validity of the third Sakharov condition. 

In [40] it was proposed that the baryon asymmetry can be generated through CP- violating 
sterile neutrino oscillations. For small Majorana masses the total lepton number of the 
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system, defined as the lepton number of active neutrinos plus the total helicity of sterile 
neutrinos, is conserved and equal to zero during the Universe's evolution. However, because 
of oscillations the lepton number of active neutrinos becomes different from zero and gets 
transferred to the baryon number due to rapid sphalcron transitions. Roughly speaking, the 
resulting baryon asymmetry is equal to the lepton asymmetry at the sphaleron freeze-out. 

The kinetics of sterile neutrino oscillations and of the transfers of leptonic number be- 
tween active and sterile neutrino sectors has been worked out in [2] . The effects to be taken 
into account include oscillations, creation and destruction of sterile and active neutrinos, 
coherence in sterile neutrino sector and its lost due to interaction with the medium, dynam- 
ical asymmetries in active neutrinos and charged leptons. For masses of sterile neutrinos 
exceeding ~ 20 GeV the mechanism does not work as the sterile neutrinos equilibrate. The 
temperature of baryogenesis is right above the electroweak scale. 

In [2] it was shown that the i/MSM can provide simultaneous solution to the problem of 
neutrino oscillations, dark matter and baryon asymmetry of the Universe. 

4- Inflation. In [3] it was proposed that the z/MSM may be extended by a light infiaton 
in order to accommodate inflation. To reduce the number of parameters and to have a 
common source for the Higgs and sterile neutrino masses the inflaton-i/MSM couplings can 
be taken to be scale invariant at the classical level and the Higgs mass parameter can be set 
to zero. The mass of the inflaton can be as small as few hundreds MeV, and the couphng 
of the lightest sterile neutrino to it may serve as an efficient mechanism for the dark matter 
production. 

5. Fine-tunings in the vMSM. The phenomenological success of the i/MSM requires 
a number of fine tunings. In particular, one of the singlet fermion masses should be in 
the 0(10) keV region to provide a candidate for the dark-matter particle, while two other 
masses must be much larger but almost degenerate [2, 40] to enhance the CP- violating 
effects in the sterile neutrino oscillations leading to the baryon asymmetry. In addition, 
the Yukawa coupling of the dark matter sterile neutrino must be much smaller than the 
Yukawa couplings of the heavier singlet fermions, to satisfy cosmological and astrophysical 
constrains [1]. These fine-tunings arc "natural" in a sense that they arc stable against 
radiative corrections. Moreover, in [41] was shown that a specific mass-coupling pattern for 
the singlet fermions, described above, can be a consequence of a lepton number symmetry, 
slightly broken by the Majorana mass terms and Yukawa coupling constants. At the same 
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time not all 18 new parameters are fixed: the allowed region in parameter space is quite 
large to yield variety of signatures to be tested with different experiments and methods. 

To summarize, none of the experimental facts, which are sometimes invoked as the argu- 
ments for the existence of the large ~ 10^° — 10^^ GeV intermediate energy scale between the 
PF-boson mass and the Planck mass, really requires it. The smallness of the active neutrino 
masses may find its explanation in small Yukawa couplings rather than in large energy scale. 
The dark matter particle, associated usually with some stable SUSY partner of the mass 
O(IOO) GeV or with an axion, can well be a much fighter sterile neutrino, practically stable 
on the cosmological scales. The thermal leptogenesis [42] , working well only at large masses 
of Majorana fermions, can be replaced by the baryogenesis through light singlet fermion 
oscillations. The inflation can be associated with the light inflaton field rather than with 
that with the mass ~ 10^^ GeV, with the perturbation power spectrum coming from inflaton 
self-coupling rather than from its mass. 

Putting all the physics beyond the Standard Model below the electroweak scale is not 
harmless, as it can be confronted with experiment at low energies (see e.g. [43] for a dis- 
cussion of neutrinoless double beta-decay in the framework of the z/MSM). The aim of this 
paper is to analyse the possibilities to search for singlet fermions responsible for baryon 
asymmetry of the Universe in the i/MSM. Finding these particles and studying their proper- 
ties in detail (in particular, CP-violating amplitudes) would allow to compute the sign and 
the magnitude of the baryon asymmetry of the Universe theoretically along the lines of [2] 
and confront this prediction with observations. The existence of the U(l) lepton symmetry 
provides an argument in favour of (9(1) GeV mass of these singlet Icptons [41]. In addition, 
the structure of their couplings to the particles of the SM is almost fixed by the data on 
neutrino oscillations. It is interesting to know, therefore, what would be the experimental 
signatures of the neutral singlet fermions in this mass range and in what kind of experiments 
they could be found. To answer this question, in this paper we will consider the variant of 
the i/MSM without addition of the inflaton; we will discuss what kind of differences one can 
expect if the light inflaton is included elsewhere. 

Naturally, several distinct strategics can be used for the experimental search of these 
particles. The first one is related to their production. The singlet fermions participate in 
all reactions the ordinary neutrinos do with a probability suppressed roughly by a factor 
(Md/MmY, where Md and Mm are the Dirac and Majorana masses correspondingly. Since 
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they are massive, the kinematics of, say, two body decays — > /jl^N, — > e^N or three- 
body decays Kl^s tt"^ + + N changes when N is replaced by an ordinary neutrino. 
Therefore, the study of kinematics of rare meson decays can constrain the strength of the 
coupUng of heavy leptons. This strategy has been used in a number of experiments for 
the search of neutral leptons in the past [44, 45], where the spectra of electrons or muons 
originating in decays of tt- and X-mesons have been studied. The second strategy is to 
look for the decays of neutral leptons inside a detector [46-49] ("nothing" leptons and 
hadrons). Finally, these two strategies can be unified, so that the production and the decay 
occurs inside the same detector [50] . 

Clearly, to find the best way to search for neutral leptons, their decay modes have to 
be identified and branching ratios must be estimated. A lot of work in this direction has 
been already done in Refs. [51-54] for the general case; we add new general results for three 
body meson decays. To analyze the corresponding quantities in the i/MSM we will constrain 
ourselves by the singlet fermion masses below the mass of the beauty mesons, Mn < 5 GeV, 
considering this mass range as the most plausible because of the reasons presented above. 
We will use the specific i/MSM predictions for the branching ratios. 

We arrived at the following conclusions. 

(i) The singlet fermions with the masses smaller than M^r are already disfavoured on the 
basis of existing experimental data of [46, 47] and from the requirement that these particles 
do not spoil the predictions of the Big Bang Nucleosynthesis (BBN) [55, 56] (s.f. [57]). 

(ii) The mass interval M-^ < Mjy < Mk is perfectly allowed from the cosmological and 
experimental points of view. Moreover, it is not excluded that further constraints on the 
couplings of singlet fermions can be derived from the reanalysis of the already existing hut 
never considered from this point of view experimental data of KLOE collaboration and of 
the E949 experiment^ In addition, the NA48/3 (P326) experiment at CERN would allow 
to find or to exclude completely singlet fermions with the mass below that of the kaon^. 
The search for the missing energy signal, specific for the experiments mentioned above, can 
be complemented by the search of decays of neutral fermions, as was done in CERN PS191 

^ We thank Gino Isidori and Yury Kudenko for discussion of these points. 
^ We thank Augusto Ceccucci for discussion of this point. 
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experiment [46, 47]. To this end quite a number of already existing or planned neutrino 
facilities (related, e.g. to CNGS, MiniBoone, MINOS or T2K), complemented by a near 
(dedicated) detector (like the one of CERN PS191) can be used^. At the same time, the 
existing setups of the MiniBooNE or MINOS experiments would unlikely allow to probe the 
cosmologically interesting parameter space of the z/MSM for Mjv < 450 MeV, where strong 
bounds on the parameters coming from CERN PS191 experiment already exist. However, 
MiniBooNE and MINOS can possibly improve the existing limits or find neutral fermions 
in the mass region 450 MeV < Mjv < M^-, where current bounds are weak (s.f. [57]). The 
record intensity of the neutrino beam at CNGS and T2K experiment are quite promising for 
heavy neutrino searches and calls for a detailed study of the possibility of neutral fermions 
detection at (possible) near detectors. 

(iii) For Mk < < Mu the search for the missing energy signal, potentially possible 
at beauty, charm and r factories, is unlikely to gain the necessary statistics and is very 
difficult if not impossible at hadronic machines hke LHC^. So, the search for decays of 
neutral fermions is the most effective opportunity. In short, an intensive beam of protons 
hitting the fixed target, creates, depending on its energy, pions, strange, charmed and beauty 
mesons that decay and produce heavy neutral leptons. A part of these Icptons then decay 
inside a detector, situated some distance away from the collision point. The dedicated 
experiments on the basis of the proton beams NuMI or NuTeV at FNAL, CNGS at CERN, 
or JPARC can touch a very interesting parameter range for Mat < 1.8 GeV. 

(iv) Going above D-meson but still below S-meson thresholds is very hard if not impos- 
sible with present or planned proton machines or B-factories. To enter into cosmologically 
interesting parameter space would require the increase of the present intensity of, say, CNGS 
beam by two orders of magnitude or producing and studying the kinematics of more than 
10^° B-mesons. 

The paper is organized as follows. In Section II we discuss the relevant part of the i/MSM 
Lagrangian and specify the predictions for the couplings of these particles coming from the 
data on neutrino oscillations and cosmological considerations. In Section III we analyze 
the present experimental and cosmological limits on the properties of these particles. In 

^ We thank Francois Vannucci for discussion of this point. 
* We thank Tasuya Nakada for discussion of this point. 
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Section IV we analyze the decay modes of singlet fermions. In Section V we consider the 
production of heavy neutral leptons in decays of K-, D- and S-mesons and of r-lepton. In 
Section 6 we analyze the possibilities of their detection in existing and future experiments. 
We conclude in Section 7. 

II. THE LAGRANGIAN AND PARAMETERS OF THE i/MSM 

For our aim it is more convenient to use the Lagrangian of the i/MSM^ in parameterization 
ofRef. [41]: 

= CusM + Niid^^^Ni - F^i L^Ni^ - MN^ ^ Ni Nj + h.c. , (1) 

where Nj are the right-handed singlet leptons (we will keep the notation without tilde for 
mass eigenstates), l>i = ^ij^*]^ ^ ^'^'^ = e,/x,T) are the Higgs and lepton doublets, 

respectively, F is a matrix of Yukawa coupling constants, M is the common mass of two 
heavy neutral fermions, AMjj are related to the mass of the lightest sterile neutrino A^i 
responsible for dark matter and produce the small splitting of the masses of N2 and A^s, 
AMjj <^ M. The Yukawa coupling constants of the dark matter neutrino |-Fqi| < lO^"*^^ are 
strongly bounded by cosmological considerations [1] and by the X-ray observations [5] and 
can be safely neglected for the present discussion and the sterile neutrino A^^i field can be 
omitted from the Lagrangian. 

In the limit AM/j — > 0, — > the Lagrangian (1) has a global U(l) lepton symmetry 
[41] . In this paper we will assume that the breaking of this symmetry is small not only in the 

^ Of course, this Lagrangian is not new and is usually used for the explanation of the small values of neutrino 
masses via the see-saw mechanism [58] . The sec-saw scenario assumes that the Yukawa coupling constants 
of the singlet fermions arc of the order of the similar couplings of the charged leptons or quarks and that 
the Majorana masses of singlet fermions are of the order of the Grand Unified scale. The theory with this 
choice of parameters can also explain the baryon asymmetry of the Universe but does not give a candidate 
for a dark matter particle. Another suggestion is to fix the Majorana masses of sterile neutrinos in 1 — 10 
eV energy range (eV see-saw) [59] to accommodate the LSND anomaly [60]. This type of theory, however, 
cannot explain dark matter and baryon asymmetry of the universe. Also, the MiniBooNE experiment [61] 
did not confirm the LSND result. The z/MSM paradigm is to determine the Lagrangian parameters from 
available observations, i.e. from requirement that it should explain neutrino oscillations, dark matter and 
baryon asymmetry of the universe in a unified way. This leads to the singlet fermion Majorana masses 
smaller than the electroweak scale, in the contrast with the see-saw choice of [58], but much larger than 
few eV, as in the eV see-saw of [59]. 
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mass sector (which is required for baryogenesis and explanation of dark matter) , but also in 
the Yukawa sector, |Fq3| <^ 1-^/92 1- For the case when iF^s] ~ |F^2| our general conclusions 
remain the same, but the branching ratios for different reactions can change. In this work 
we also neglect all CP-violating effects, which go away if the lepton number symmetry is 
exact. 

To characterize the measure of the U{1)l symmetry breaking, we introduce a small param- 
eter e = F3/F < 1, where = [F^fFja, and F2 = F. As was shown in [41] , there is a lower 
bound on e coming from the baryon asymmetry of the Universe, e > 0.6 • k • 10~^(M/GeV), 
where k = 1(2) for the case of normal (inverted) hierarchy in active neutrino sector. 

The mass eigenstates (iV2,3 without tilde) are related to N2^3 by the unitary transforma- 
tion, 

N^UrN, (2) 
where the 2x2 matrix Ur has the form 



V2 



(3) 



where the phases 0^ can be expressed through the elements of AMjj, the explicit form of 
which is irrelevant for us. 

As a result, for e <S 1 the interaction of the mass eigenstates N2 and A^3 has a particular 
simple form, 

1 - -Mo - Mq - 

Lat ~ --^fMN2 + N^)^ - -fN2'N2 - -fN^'^N^ + h.c. , (4) 

where = \Fa.2\- The masses M2 and M3 must be almost the same (baryogenesis con- 
straint), AM^ = |M| — M|| < 10~^M^ [2, 40, 41]. The baryon asymmetry generation occurs 
most effectively if AM^ ~ (2 keV)^, but smaller and larger degeneracy works well also. 

The fact that two heavy fermions are almost degenerate in mass may be important 
for analysis of the experimental constraints. In decays of different mesons or r-lepton a 
coherent combination (A^2 + ^3) will be created, while in a detector of size / situated on the 
distance L from the creation point an admixture of the (A^2 — -^3) state with the probability 
(in the relativistic limit) P ~ sin^ (j) will appear {E is the energy of the neutral fermion, 
(f) — AM^L/(4£')). For ^ 1 coherence effects are not essential and the description of 
the process in terms of A^2 and A^3 is completely adequate, while if ~ 1 the coherence 



effects are important, and order e terms describing tlie interactions of (A^2 — N^) witfi tfie 
particles of the SM must be included. Numerically, if AM^ >(2 keV)^, Z ~ 10 m, and 
E < 100 GeV, then (f)l/L> 10^, and N2 oscillations can be safely neglected. Only this 

case will be considered in what follows. 

As it was demonstrated in [41], the coupling constants can be expressed through the 
elements of the active neutrino mass matrix M^. To present the corresponding relations, we 
parameterize Mi, following Ref. [62]: 

M^^V* ■ diag(mi, m2e^'^' , m^e^'^^) ■ , (5) 

with V = -R(^23)diag(l, e^''^, l)R(^^i3)R(^^i2) the active neutrino mixing matrix [63], and 
choose for normal hierarchy mi < 1112 < and for inverted hierarchy 1713 < mi < m2- AH 
active neutrino masses are taken to be positive. As was shown in [1, 4-6], the one of the active 



neutrino masses must be much smaller than the solar mass difference, mgoi = s/ ^^soi — O-^l 
eV, so that other active neutrino masses are simply equal to matm = •>/ Am^tm — 0.05 eV and 
to msoi for the case of normal hierarchy and to matm with a mass splitting 5m — ml^j2mg,tm 
for the case of inverted hierarchy. 
The coupling F is given by [41]: 



2 m^tmM 

where v — 174 GeV is the vev of the Higgs field and k, ~ 1(2) for the case of normal 
(inverted) hierarchy. 

The ratios of the Yukawa couplings can be expressed through the elements of the 
active neutrino mixing matrix [41]. A simple expression can be derived for the case ^13 = 
0, O2Z — 'k/A, which is in agreement with the experimental data. For normal hierarchy there 
are possibilities: 

fe ■■ fl ■■ fr « —sm'e,2\l ±xf:hl-x'\':l\l±x\\ (7) 
m3 z z 



where x = ie^^^'^ ^'^ ^^^^ /— cos ^12, and all combinations of signs are admitted. For a nu- 
merical estimate one can take [62] sin^ ^12 — 0.3, leading to a; ~ 0.35ie*('^^~''^~'^^^ and to 
fe/if^ + /^) ~ 0.05. In other words, the coupling of the singlet fermion to the leptons of 
the first generation is suppressed, whereas the couplings to the second and third generations 
are close to each other. 



10 



For the case of inverted hierarchy two out of four solutions are almost degenerate and 
one has [41] : 

fe fl fr ~ (8) 

where p — ± sin (^1 sin(2^i2). Taking the same value of 612 as before, we arrive at fe/{f^ + 
f^) ~ (0.04 — 25), depending on the value of unknown CP- violating phase Si. The couplings 
of A^2,3 to II and r generations are almost identical, but the coupling to electron and its 
neutrino can be enhanced or suppressed considerably. The corrections to relations (7,8) are 
of the order of 0{e) and for e ~ 1 the ratios of the coupling constant can be quite different 
from those in eqns. (7,8). 

The relations (6,7,8) form a basis for our analysis of experimental signatures of heavy 
neutral leptons. In most of the works the strength of the coupling of a neutral lepton X to 
charged or neutral currents of flavour a is characterized by quantities Uax and Vax- In the 
case of the i/MSM there are two neutral leptons with almost identical couphngs (if e <S 1), 
so that 

|C/al|-|Kl| = |^a2| = |K2| = |^a| • (9) 

The overall strength of the coupling is given by 

a 

whereas the relations between different flavours follow from (7,8). 

As it was found in [41] (see also [2, 40]), for successful baryogenesis the constant F must 
be small enough, F< 1.2 x 10~^, otherwise N2 and A^3 come to thermal equilibrium above 
the electroweak scale and the baryon asymmetry is erased. This leads to the upper bound 

< 2k X 10-' i^-j^j . (11) 

It is the smallness of the required strength of coupling which makes the search for neutral 
leptons of the vMSM be a very challenging problem, especially for large M. 

In the framework of the z/MSM, a lower bound on U can be derived as well. The maximal 
value of the parameter e, characterizing the breaking of the U(l) leptonic symmetry is e = 1. 
This results in 

> l-3« X 10-- . (12) 
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Further cosmological constraints on the couphngs of heavy sterile neutrinos are coming 
from BBN. The cosmological production rate of these particles peaks roughly at the tem- 
perature [55] Tpeak ~ 10(M/GeV)^/^ GeV and for C/^ > 2 x IQ-^^ ^GeV/M) they were 
in thermal equilibrium in some region of temperatures around Tp^ak- This is always true, 
since in the z/MSM the constraint (12) is required to be valid. We will see below that the 
BBN constraints are in fact stronger than those of (12) for relatively small fermion masses 
Mn < 1 GeV. On the basis of inequalities (11,12) and hmits from BBN, the i/MSM can be 
probed (either confirmed or ruled out) in particle physics experiments. 

The relations (7,8) still allow a lot of freedom in relations between Yukawa couplings to 
different leptonic flavours, since the Majorana CP-violating phases in the active neutrino 
mass matrix are not known. Therefore, to present quantitative predictions we will consider 
three sets of Yukawa couplings corresponding to three "extreme hierarchies" , when value 
of Yukawa constants fa, fp are taken to be as small as possible compared to another one 
f^, a ^ (5 ^ which thus mostly determines the overall strength of mixing C/^. In what 
follows we will refer to these sets as benchmark models I, II and III with ratios of couphng 
constants which can be read off from eqs. (7), (8): 



model I : 


fl 


■ f 


■■fr 


Ri 52 : 1 : 


: 1 , K= 1 , 


model II : 


fl 




■■fr 


^ 1 : 16 : 


:3.8, K = 2, 


model III : 




■ f : 


■■fr 


fa 0.061 : 


1 : 4.3 , K = 2 



Let us explain how these numbers were obtained. For the model I we simply increase 
in a maximal way the value of the coupling constant to electron, choosing the appropriate 
combination of signs in eq. (8). In case of model III the couphng of N to the third generation 
of leptons is stronger than to the others. This could only happen if the hierarchy of active 
neutrino masses is normal, see eq. (7). Choosing real and positive x one can see that the 
maximum value of the ratio l/r/Z/xP is given by 

As reference point we choose the central values of parameters of neutrino mixing (see, e.g. 
[62]), that gives x ^ 0.35. This means that the ratio (13) can be as large as 4.3 (varying 
the parameters of the active neutrino mixing matrix within their error bars one arrives at a 
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bit larger number). By the same type of reasoning the maximal values of the ratio |/r//e|^ 
is given by 

sin'^i2-( M^^l I -71. (14) 




I/el' \2m3 

Similar considerations provide values of Yukawa couplings in model II. 

These benchmark models are choosen to show the variety of quantitative predictions 
within originally 18-dimcnsional parameter space of z/MSM, constrainted already by cos- 
mology, astrophysics, and observations of neutrino oscillatuions. For a given process, they 
should be confined between numbers given for benchmark models for e <^ 1. A special study 
should be undertaken to outline the actual range of i/MSM predictions in case of e ~ 1, when 
relations (7) and (8) become invalid. 



III. LABORATORY AND BBN CONSTRAINTS ON THE PROPERTIES OF 
HEAVY LEPTONS 

The aim of this section is to discuss whether the past experiments devoted to the search 
for neutral leptons have entered into cosmologically interesting parameter range defined by 
eqns. (11,12). In addition, we will consider the Big Bang Nucleosynthesis constraints on the 
properties of heavy leptons in the i^MSM. 

The analysis of the published works of different collaborations reveals that for the mass 
of the neutral lepton M > 450 MeV none of the past or existing experiments enter into 
interesting for i/MSM region defined by eq. (11). The NuTeV upper limit on the mixing 
is at most 10"^ in the region M ~ 2 GeV [48], whereas the NOMAD [49] and L3 LEP 
experiment [50] give much weaker constraints. Note that the eqns. (11,12) give at M = 2 
GeV: 6 • 10-^' < U^/k < 5 • lO'^. 

The best constraints in the small mass region, M < 450 MeV are coming from the CERN 
PS191 experiment [46, 47], giving*^ roughly |[/e,Mp ^ 10"^ the region 250 MeV < M < 450 

^ The most recent published results of CERN SPS experiment [47] contain the exclusion plots up to 400 
MeV. In a previous publication, [46], the limit on U^, though not as strong as in [47], was presented up 
to 450 MeV. We became aware of PhD Thesis of J.-M. Levy [64] (we thank F. Vannucci for providing us 

a copy of this manuscript) which contains the experimental exclusion plots for and \UeUf^\ up to 450 
McV. Wc use these unpublished results in our work. If the results of [64] arc ignored, our plots should 
be modified accordingly in the region 400 MeV < Mjv < 450 MeV, and phenomenologically viable region 



13 



MeV (the NuTeV limit in this mass range is some two orders of magnitude weaker) . These 
numbers are already in the region (11) and thus provide non-trivial limits on the parameters 
of the i/MSM. Moreover, as it will be seen immediately, the considerations coming from 
BBN allow to establish a number of lower bounds on the couplings of neutral leptons which 
decrease considerably the admitted window for the couplings and masses of the neutral 
leptons. 

The successful predictions of the BBN are not spoiled provided the life-time of sterile 
neutrinos is short enough. Then neutrinos decay before the onset of the BBN and the 
products of their decays thermalize. This question has been studied in [55] and we will use 
the results of their general analysis for the case of Models 1-111 described in Section 11. 

First, we note that [55] considered the case of one sterile neutrino of Dirac type, whereas 
we have two Majorana sterile neutrinos^. This means that we have exactly the same number 
of degrees of freedom and that the constraints of [55] , expressed in terms of lifetime of sterile 
neutrino are applicable to our case. 

Ref. [55] studied in detail only the mass range 10 MeV < Mjv < 140 MeV, for higher 
masses these authors argued that the life-time tm of the heavy lepton must be smaller than 
0.1 s to definitely avoid any situation when heavy lepton decay products could change the 
standard BBN pattern of light element abundances. We note in passing that it would be 
extremely interesting to repeat the computation of [55] for Mjv > 140 MeV in order to have 
a robust BBN constraints in this mass range; meanwhile we will just require (conservatively) 
that Tiv < 0.1 s for neutral fermions heavier than 7r-meson. 

For the masses in the interval 10 MeV < M < 140 MeV the constraint on the mixing 
angle, based on a fit to numerical BBN computations [55], reads 



with si,e = 140.4, si,^ = si,r = 568.4, S2,e = -1-05 • 10-^ S2,^ = S2,r = "5.17 • 10-6, 
= —3.070 and = a-r = —3.549 (we took a conservative bound equivalent to adding 
one extra neutrino species, as explained in [55]); the limits (15) are valid in the models where 
sterile neutrino mix predominantly with only one active flavor. Here we took into account 



The concentration of the dark matter sterile neutrinos is weU below the equilibrium one so that its existence 
may be safely neglected at this time. 




(15) 



expands. 
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that in Ref. [55] neutrinos of Dirac type have been considered, while we discuss neutrino 
of Majorana type, hence the total width contains an extra factor 2 in comparison with the 
Dirac case and the constraint of is in fact 2 times weaker than that of [55]. The limits 
(15) can be converted into limits on the mixing f/^ for the models 1-111. 

To consider higher masses we computed the life-time of heavy leptons (the details of 
computation can be found in Section IV) and required that it exceeds 0.1 s, to make a 
conservative exclusion plot. The most important decay channels for Mjv < Mk are the 
two-body semileptonic ones N — > n^u, N — > n^e"^, N — > n^n^. 

For various patterns of neutrino mixing we present the experimental and BBN constraints 
in Fig. 1. Note that in extracting the limits on mixing from [46, 47] (this experiment 
presented 90% confidence level exclusion plot) we also take into account that there arc two 
degenerate neutrinos in the z/MSM, and that the constraints in [46, 47] are given for Dirac 
type sterile neutrinos. For the same value of the mixing angles, the same number of sterile 
neutrino helicity states are created in both Dirac and Majorana cases, but in the former 
case only half of states contribute to each decay channel. Hence, the constraints on |C/eP, 
|f/e||?7^| and |f/^P are in fact by a factor 2 stronger, since the number of decay events is 
proportional to \U\'^. 

One can see that depending on the type of the neutrino mass hierarchy and specific 
branching ratios in the benchmark models I-III the phenomenologically allowed region of 
parameter space can be reduced or enlarged. Moreover, the masses below the tt meson mass 
are excluded in most cases^ but still there are models where small regions of the parameter- 
space above the pion mass are perfectly allowed^. We would also like to stress that the 
branching ratios for e ~ 1 can be quite different from (7,8) leading to extra uncertainties. 

Above pion mass, the BBN limits are down to two order of magnitude below the direct 
limits form CERN PS191 experiment, thus one-two orders of magnitude improvement is 
required to either confirm or disprove the i/MSM with sterile neutrinos lighter than 450 MeV. 
For the three benchmark models we transfered these limits to the upper limits on overall 

^ For the i/MSM with hght inflaton BBN bounds are weaker and masses below pion are certainly allowed 
[41]. 

^ Note that our exclusion plot is different from that of Rcf . [57] , where the coupling of sterile neutrino to r 
generation was not considered. Moreover, eq. (3.1) of this paper contains a factor 4 error. In addition, 
the formula (21) of [56] for the probability of N ^ w^v decay is not correct, see discussion in Section IV. 
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FIG. 1: Limits on |C/eP) |f^e||t^^| and \U^\'^ for three benchmark models (I-III from left to right) 
from BBN (lower bound) and from direct searches in the CERN PS191 experiment (upper bound). 
Blank regions are phenomenologically allowed. 

mixing f/^ and neutrino lifetime and plotted them in Fig. 4. 

The improvement required to test the i/MSM with sterile neutrinos lighter than 450 MeV 
can be done with either new kaon experiments, such as one planned in JPARC, or special 
analysis of the available data on kaon decays collected in Brookhaven and Frascati. In 
particular, E787/E949 Collaboration reported limit on —>■ tt+X decay with X being 
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hypothetical long-hved neutral particle [65] . With statistics of thousand of billions charged 
kaons, available in this experiment, one can expect to either prove or completely rule out 
vMSM with sterile neutrinos lighter than 450 MeV. The same conclusion is true for the third 
stage of CERN NA48 experiment. 

In the next two Sections wc discuss the decays and production of neutral fermions for 
a mass range up to 5 GeV, to understand the requirements to possible future experiments 
that could allow to enter into interesting parameter space for neutral fermion masses above 
400 MeV. 

IV. DECAYS OF HEAVY NEUTRAL LEPTONS 

Heavy neutral leptons wc consider {Mn ^ 10 MeV) are unstable, since decay channels to 
light active leptons, N — > I'a^a^is, ^ ~^ e~^e~i'a are open; the modes like A^2,3 ^ A^i + . . . 
are strongly suppressed. Hereafter charge conjugated modes are also accounted resulting in 
double rates for Majorana neutrinos as compared to Dirac case. For heavier leptons more 
decay modes are relevant, 

N — > fxeu, n^u, ne, ji^ jjTv^ nfj,, Ke, Kjj, , 771/, pi/, . . . 

Decays of sterile neutrinos have been exhaustedly studied in literature. For convenience 
we present explicit formulae for relevant decay rates in Appendix A. Most of them (but not 
all) can be be obtained straightforwardly by making use of the formulae for Dirac neutrinos 
presented in Ref. [54], which we found to be correct. 

Neutrino decays branching ratios for benchmark models I-III and M^- < Mj^ < 2 GeV 
are plotted in Figs. 2, 3. For heavier neutrino many-hadron final states become important, 
and one can use spectator quarks to calculate the corresponding branching ratios. Below 
2 GeV the contribution of these modes to total neutrino width is less than 10%. Neutrino 
hfetime is constrained by limits (11), (12) on overall strength of mixing. The results for 
models I, II and HI are presented in Fig. 4a: in phenomenologically viable models neutrino 
lifetime is confined by corresponding solid (upper limits) and dashed (lower limits) lines. 
The horizontal solid line indicates the order-of-magnitudc upper limit on neutrino life time. 
Tat < 0.1 c, which guarantees that the results of standard BBN remain intact [55] for 
Mjv ^ 140 MeV. In a given model the range of neutrino mass, where the corresponding solid 
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FIG. 2: Branching ratios of neutrino two-body decays Ni — > XY as functions of neutrino mass 
Mjv for models with the same hierarchy in mixing as in models: a) I, b) II, c) III; different lines 
correspond to different modes: irv (solid black), 7re (short-dashed black), tt/x (long-dashed black), 
Ke (solid light gray), 'qv (solid dark gray), rj'v (short-dashed dark gray), Kjjl (dashed light gray), 
pu (solid gray), pe (short-dashed gray), pp (long-dashed gray). 

linc(s) is(arc) above the corresponding dashed one(s) is disfavoured. 

These limits imply limits on overall mixing U"^ plotted in Fig. 4b: in phenomenologically 
viable models mixing U"^ is confined by corresponding solid and dashed lines. One can see 
that the constraint from BBN is stronger than the see-saw constraint (12) for M < 1 GeV. 
However, it is worth noting that the limit tjv < 0.1 s may happen to be too conservative and 
can presumably be relaxed to some extent provided careful study of processes in primordial 
plasma in BBN epoch. In what follows, for the three benchmark models we give upper and 
lower limits on various neutrino rates. For a given neutrino mass these limits are saturated 
respectively by the tightest among upper limits and tightest among lower limits on neutrino 
mixing, presented in Fig. 4b. Only these tightest limits are used below. 

Note in passing that as we already mentioned the i/MSM predictions beyond benchmark 
models could deviate to some extent from a naive interplay between benchmark numbers. 
At the same time for any set of parameters the presence of both upper and lower bounds on 
neutrino rates is a general feature of i/MSM, which allows it to be falsified. 
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FIG. 3: Branching ratios of neutrino three-body decays Nj — > ABC as functions of neutrino mass 
Mjv for models with the same hierarchy in mixing as in models: a) I, b) II, c) III; different lines 
correspond to different modes: uui/ (sum over all invisible modes, solid black), z/e+e~ (solid gray), 
ue^ (sum over two modes, long-dashed gray), vjjl^ix" (short-dashed gray). 

V. PRODUCTION OF HEAVY NEUTRAL LEPTONS 

In high-energy experiments the most powerful sources of heavy neutral leptons are the 
kinematically allowed weak decays of mesons (and baryons) created in beam-beam and 
beam-target collisions. Obviously, the relevant hadrons are those which are stable with 
respect to strong and electromagnetic decays. 

The spectrum of outgoing heavy neutral leptons in a given experiment is determined 
mostly by the spectrum of produced hadrons H subsequently decaying into heavy leptons. 
Since relevant hadrons contain one heavy quark Q, differential cross section of their direct 
production dcr^'' can be estimated by use of the factorization theorem 

— — ^;2— =/ dz-5(pQ-zpH)-DH,Q(z)-— — |^ , (16) 
dpH,LdpHT Jo dpQ^LdpQT 

where dag^ is differential cross section of direct Q-quark production^*^, ph,l, Ph,t and pq,l, 
Pq^t are longitudinal and transverse spatial momenta of hadron H and heavy quark Q, 

We assume non-polarized beam(s) and target and hence axial symmetry. 
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FIG. 4: a) Upper (solid lines) and lower limits (dashed lines) on neutrino lifetime in models I 
(black), II (dark gray) and III (light gray); a horizontal thick solid black line indicates the upper 
limit from BBN, rjv < 0.1 c, suggested in Ref.[55], thin solid lines are limits from eq. (12), thick 
dashed lines refer to eq. (11), thin dashed lines correspond to limits from direct searches for sterile 
neutrinos discussed in Section III; charge-conjugated modes arc accounted, b) Lower (solid lines) 
and upper limits (dashed lines) on overall mixing C/^ in models I (black) , II (dark gray) and III (light 
gray). Thin solid lines and thick dashed lines depict limits from eqs. (12) and (11), respectively. 
Thick solid lines indicate lower limits from order-of-magnitude BBN bound on neutrino lifetime, 
ttv < 0.1 c for Mn > 140 MeV, thin dashed lines refer to limits from direct searches for sterile 
neutrinos discussed in Section III. 

respectively; zpn is a part of hadron momentum carried by heavy quark and a fragmenta- 
tion function Dh,q{z) describes the details of hadronization. The differential cross section 
entering the integrand in eq.(16) can be calculated within perturbative QCD, while func- 
tion Dh,q{z) comprises non-perturbative information. There are several approximations to 
Dh,q{z) in literature, e.g. commonly used in high energy physics generator PYTHIA adopts 
modified Lund fragmentation function [66] 



D(z) OC — — , 



with default parameters a — 0.3 and b — 0.58 GeV 
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The rate of hadron production depends on the intensity of coUisions. The distribution of 
total number of directly produced hadrons dNf}'^ reads 



dpH,LdpH,T dpH,LdpH^T 

where Cacc is an integrated luminosity of a given experiment and we neglect tiny imprints of 
real bunch structure on outgoing hadronic spectra. Note that we are interested in hadrons 
stable with respect to strong and electromagnetic decays, thus apart of direct production 
they emerge due to strong and electromagnetic decays of other hadrons, which give indirect 
contribution dN^^'^. The distribution of the total number of produced hadrons dN^ is a sum 
of both contributions, 

dNii _ dN'll'' ^ dNf^'^ 



dpH,LdpHT dpH,LdpjjT dpH,LdpjjT 

Produced hadrons stable with respect to strong and electromagnetic decays travel dis- 
tances of about • th • jH {Ph, th and are speed, lifetime and boost factor of a 
given hadron) and then decay weakly, producing some amount of heavy neutral leptons. In 
the hadron rest frame the spatial momentum of heavy lepton p^ can be correlated with 
the hadron total spin. Consequently, in the laboratory frame there can be additional to 
Lorenz boost contribution to correlations between pj^ and pu- This contribution is smearing 
with growth of statistics and can be also neglected if typical 7-factor of hadrons is large, 
7h = Eh/Mh 3> 1. Hence, in the laboratory frame, the distribution of heavy leptons over 
spatial momentum is given by 



dNy _ r dBn (H ^N + ...) 

dpN,Ldpjf^T H " ^ 



X 



J d^n^ ■ S ^Pn - Ph - • \J Ej^ - Ml 



dN 



(17) 



11 



dpH,LdpjjT 



where we integrate over unit sphere boosted to laboratory frame and sum up all contributions 
from all relevant hadrons; dBn {H ^ N + . . .) is a differential inclusive branching ratio of 
hadron H into heavy neutrino. These branching ratios can be straightforwardly obtained 
for each hadron with help of the standard technique used to calculate weak decays in the 
framework of the SM. Indeed, in both models (MSM and i/MSM) neutrinos are produced 
mostly via virtual VF-boson (charged current): the only difference is that in i/MSM neutrinos 
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are massive. For heavy neutrinos this results^^ in enhancement of pure leptonic decay modes 
which are strongly suppressed in the SM by charged lepton masses. 

The heavier the quark the lower its production rate; hence, a class of the lightest kinemat- 
ically allowed hadrons saturates heavy neutrino production. As we explained in Section III, 
neutrinos in phenomenologically viable i/MSM are likely to be heavier than pion. If neutrino 
Ni is lighter than kaon, the dominant source of neutrinos is decaying kaons, 

^ l^Nj , (18) 

Kl ^ TT^l^Nj . (19) 

The two-body decays (18) have been already studied in literature (see, e.g., Refs. [52, 53]). 
For convenience, the differential branching ratio is presented in Appendix B. Contribution 
of three-body decays (19) to neutrino production is suppressed by phase volume factor; as 
the largest impact they give a few per cent at Mjv ^ M^; the corresponding differential 
branching ratio is presented in Appendix B. 

In models with neutrino heavier than kaon but lighter than charmed hadrons, decays of 
those latter dominate neutrino production. The largest partial width to heavy neutrinos is 
exhibited by D^-meson which leptonic decays Dg — > laNj are not suppressed by CKM mixing 
angles as compared to similar decays of D-mesons, — > I'^Nj. Semileptonic three-body 
decay modes 

^ ri^'^l^Ni , D ^ Kl^Ni , (20) 

Ds -> (l>laNj , D ^ K*laNj (21) 

are unsuppressed by CKM-mixing as well, and are sub-dominant in general. For sufficiently 
hght neutrinos, Mk ~ < 700 MeV, £)-meson semileptonic decay modes give contri- 
bution comparable to Dg la^i at Mk ^ Mjv ^ 700 MeV, because the D-meson total 
production dominates over D^-production in hadronic collisions. Differential branching ra- 
tios of the leptonic decays and the semileptonic decays (20) of charmed mesons are provided 
by general formulae in Appendix B, where the expression of differential branching ratio to 
vector mesons V — (t),K* (21) is also presented. Both the rest of kinematically allowed 

Also, in models with heavy neutrinos values of hadronic form factors governing semileptonic width are 
changed in accordance with shift in virtuality of W-boson. 



22 



three-body decay modes and four-body decay modes, e.g. D — > KnlaNi, are strongly sup- 
pressed by either CKM mixing or phase volume factor and can be neglected. The largest 
contribution from charmed baryons comes from the decay Ac — > MaN and is negligibly small 
for heavy neutrino production. 

In models where neutrino masses are within the range 2 GeV ^ M^v ^ 5 GeV, neutri- 
nos are produced mostly in decays of beauty mesons. These are also mostly leptonic and 
semileptonic decays, which branching ratios are described by general formulae presented in 
Appendix B. As compared to £)- meson decays, S-meson decays into heavy neutrinos are 
strongly suppressed by off-diagonal entries of CKM matrix. For neutrinos hghter than about 
2.5 GeV semileptonic modes to charm mesons, e.g. B — > D^*HNi, dominate over leptonic 
mode B INj because of both larger CKM- mixing, \Vi,c\ ^ and larger values of 

hadronic form factors, Jb/Mb ^ /+,/o- Be INj is more promising, but Sc-production 
in hadron collisions is suppressed. For heavier neutrinos leptonic modes dominate. The 
baryon contribution is subdominant at any Mat. 

Note, that additional, but always subdominant, contribution to heavy neutrino produc- 
tion comes from decays of r-leptons (if kinematically allowed), which emerge as results of 
decays of Dg- and 5-mesons. 

The total number of produced heavy leptons is given by the integration of eq. (17) 
over iLy and Ej^. For order-of-magnitude estimates one can use the following simple approx- 
imation, 

^iv = J]iV^-Br(//^Ar...) , 

H 

with Nff being a total number of produced hadrons H, which in turn can be estimated as 

Nn^NQ-BriQ^H) , 

where Nq is a total number of produced heavy quarks Q and Br {Q H) isa relative weight 
of the channel Q ^ H in Q-quark hadronization. For strange meson the reasonable estimate 
is Bi:{s^ K-) = Br {s Kl). Following Ref. [67] we set Br (c D+) = 0.4 • Br (c ^ Dq) 
and assuming Br (c Dg) = Br (c ^ Ac) obtain for relevant hadrons 

Br (c ^ £>+) = 0.2 , Br (c ^ Do) = 0.5 , Br (c ^ L>^) = 0.15 . 

This is a consequence of strong overlapping between quark wave functions in the meson required to produce 
virtual T4^-boson in case of leptonic decay. 
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For beauty mesons we use [68] 

Br {b B+) = Br (6 ^ = 0.4 , Br (6 ^ Bs) = 0.1 . 

For each heavy quark Q the dominant contribution to heavy neutral lepton production 
comes from leptonic and semileptonic decays of mesons. The hmits on branching ratios for 
relevant decays are plotted in Figs. 5-15 as function of neutrino mass for three benchmark 
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FIG. 5: Branching ratios of decays K — > eNi (solid lines) and K jjbNi (dashed lines) as 
functions of heavy neutrino mass in models: a) I, b) II, c) III. In a phenomenologically viable 
model and heavy neutrino mass within ^ ^ Mk, the branching ratios are confined 
between corresponding thin and thick lines which show upper and lower limits on U"^ from Fig. 4b, 
respectively. 

models. Within i/MSM the interesting branching ratios are confined between corresponding 
thin (upper limit) and thick (lower limit) lines: inside these regions all limits on f/^ plotted 
in Fig. 4 are fulfilled, in a given model the neutrino mass region, where the corresponding 
thin line is below the corresponding thick line, is disfavoured. Rate doubling due to heavy 
neutrino degeneracy is taken into account. 
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FIG. 6: Branching ratios of decays D eNj (gray solid lines), D fiNj (gray long-dashed 
lines), Dg — eNj (black solid lines), Dg — ^Nj (black long-dashed lines) and Dg — tNj (black 
short-dashed lines) as functions of heavy neutrino mass Mjv in models: a) I, b) II, c) III. In a 
phenomenologically viable model and heavy neutrino mass within M^^ < Mjv ^ Md, the branching 
ratios are confined between corresponding thin and thick lines which show upper and lower limits 
on U'^ from Fig. 4b, respectively. 

The two-body decays can be searched for to probe i/MSM: produced charged leptons are 
monochromatic with spatial momenta 



The positions of these peaks in charged lepton spectra and their heights are correlated 
obviously for different modes and mesons. These features is a very clean signature of heavy 
leptons. From the plots in Fig. 6 one concludes that statistics of billions charmed hadrons 
is needed to probe z/MSM with neutrino of masses 0.5 GeV ^ Mn ~ 2 GeV. In models with 
lighter neutrinos kaon decays are important and required statistics is smaller. Contrary, in 
models with heavier neutrinos statistics has to be larger and it is a challenging task for future 
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FIG. 7: Branching ratios of decays B eNj (black solid lines), B fxNj (black long-dashed 
lines), B — tNj (black short-dashed lines). Be eNj (gray solid lines). Be fiNj (gray long- 
dashed lines) and Be — tNj (gray short-dashed lines) as functions of heavy neutrino mass Mjv in 
models: a) I, b) II, c) III. In a phenomenologically viable model and heavy neutrino mass within 
^ Mat ^ Mb, the branching ratios are confined between corresponding thin and thick lines 
which show upper and lower limits on C/^ from Fig. 4b, respectively. 

5-factories. Note that the set of phenomenologically interesting models where neutrinos are 
produced in kaon decays can be examined completely, as it requires billions of kaons and 
collected World statistics is much larger. 

Semileptonic decays also contribute to heavy lepton production, but spectra of outgoing 
leptons and mesons are not monoenergetic, making this process be less promising probe of 
i/MSM heavy neutrinos. 

To illustrate the relative weight of different mesons in total neutrino production we plot 
in Fig. 16 the quantity 

= J2^Q,h ' ^Q^H ^Br{Q^H)-Br{H^N...) 

H 

(where all considered above leptonic and semileptonic decays of strange, charmed and beauty 
mesons are taken into account, Q — s,c,b) within relevant ranges of neutrino masses Mn. 
With a reasonable estimate of strange, charm and beauty cross sections at large ener- 
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FIG. 8: Branching ratios of semileptonic decays K ireNj (solid lines) and K ttjiNi (dashed 
lines) as functions of heavy neutrino mass Mn in models: a) I, b) II, c) III. In a phenomenologically 
viable model and heavy neutrino mass within ~ Mn ^ Mk, the branching ratios are confined 
between corresponding thin and thick lines which show upper and lower limits on from Fig. 4b, 
respectively; form factors are taken from Refs. [68]. 
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one concludes that to produce a few neutrinos lighter than kaon, 10^-10^'' collisions is 
required, while for heavier neutrinos the statistics should be four orders of magnitude 
(0.5 GeV < < 2 GeV) or even eight orders of magnitude (2 GeV < < 4 GeV) 
larger. 

Note in passing that in our considerations baryon decays as well as decays with more 
than three particles in a final state have been neglected. These additional contributions to 
neutrino production are expected to be insignificant. 
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FIG. 9: Branching ratios of semileptonic decays D ireNj (black solid lines), D irfiNj (black 
dashed lines), D KeNj (gray solid lines) and D — > KjiNi (gray dashed lines) as functions of 
heavy neutrino mass Mjv in models: a) I, b) II, c) III. In a phenomenologically viable model and 
heavy neutrino mass within M^r Mn < Md, the branching ratios are confined between corre- 
sponding thin and thick lines which show upper and lower limits on from Fig. 4b, respectively; 
form factors are taken from Refs. [69]. 

VI. PROSPECTS FOR FUTURE EXPERIMENTS 

Generally, there are two types of processes where heavy neutrinos can be searched for: 
neutrino production hadron decays and neutrino decays into SM particles. 

In Section V we presented plots with hadron branching ratios to neutrinos in the frame- 
works of the three benchmark models. ^From these plots one can conclude that statistics 
expected at proposed Super B-factories give a chance to explore i/MSM with neutrinos hghter 
than about 1 GeV and probe some part of parameter space, if neutrino masses are in 1- 
2 GeV range. For heavier neutrinos typical branching ratios become too small, so even with 
large number of available hadrons actually small uncertainties in prediction of background 
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FIG. 10: Branching ratios of semileptonic decays Dg XlNj, X = r],rf',K (black, dark gray, 
light gray lines), I = e,iJ, (solid and dashed lines), as functions of heavy neutrino mass in 
models: a) I, b) II, c) III. In a phenomenologically viable model and heavy neutrino mass within 
< Mjv ^ Md, the branching ratios are confined between corresponding thin and thick lines 
which show upper and lower limits on C/^ from Fig. 4b, respectively; form factors are taken from 
Ref. [70]. 

can make any searches insensitive. 

For heavier neutrinos the most promising experiments are beam-target experiments with 
high intensity of a beam and high energy of incident protons. Heavy neutrinos from decays 
of numerous secondary hadrons will travel some distance and then decay into SM particles 
with branching ratios discussed in Section IV. With lifetime in the range 10~^ -j- 10~^ s 
neutrino covers a distance in exceed of one kilometer, so a detector aimed at searches for 
neutrino decay signatures should be placed at an appropriate small distance from the target 
to avoid decrease of statistics due to neutrino beam divergence. In what follows we consider 
the experimental setup with appropriately thin target, assuming that produced in beam- 
target collision hadrons decay freely without further interaction inside the target. So, this 
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Mat, GeV 



Mjv, GeV 



Mjv, GeV 



FIG. 11: Branching ratios of semileptonic decays D K*eNi (black solid lines), D K*iiNi 
(black dashed lines), Dg — (t)eNi (gray solid lines) and Dg — ^/xA?/ (gray dashed lines) as functions 
of heavy neutrino mass Mat in models: a) I, b) II, c) III. In a phenomenologically viable model and 
heavy neutrino mass within M-,^ ^ M^v ^ -/^d, the branching ratios are confined between corre- 
sponding thin and thick lines which show upper and lower limits on from Fig. 4b, respectively; 
form factors are taken from Refs. [70, 71]. 

is not a classical beam-dump setup. For classical beam-dump experiment secondary kaons 
interact in material before decay, that change their contribution to production of neutrinos 
with Miv < Mk-i which estimate requires additional study. Heavier neutrinos are produced 
mostly by D- and S-mesons, which even in beam-dump setup decay before interaction. 
Hence, for > Mk our results obtained below are valid for beam-dump experiment as 
well. 

The total number of neutrinos produced by Npor incident upon a target protons with 
energy E is given by 
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J J. X.J \J \J . -L -l-.J - 

Mtv, GeV Mn, GeV Mjv, GeV 



FIG. 12: Branching ratios of decays r — ttA/" (black solid lines), r — ifiV (black dashed lines), 
r — pN (dark gray solid lines) , r — ^ z/eA^/ (sum over all active neutrino species, dark gray dashed 
lines), T ffiNi (sum over all active neutrino species, light gray solid lines) as functions of heavy 
neutrino mass Mjv in models: a) I, b) II, c) III. In a phenomenologically viable model and heavy 
neutrino mass within M^r ^ Mt, the branching ratios are confined between corresponding 

thin and thick lines which show upper and lower limits on from Fig. 4b, respectively. 

where A refers to the target material and Mpp{E) is a total multiplicity (average number 
of secondary particles in proton-proton collision). Here we suppose that all beam protons 
interact once in the target; the account of finite thickness of the target is straightforward 
and results in effective decrease in NpoT- Assuming as a reasonable approximation at large 
E 

we arrived at 

Nr,{E) = J2 ^Q- ^q(^) ■ NpoTiE) ■ Mpp{E) . 

Q=s,c,b 

Below we present the numerical estimates for four high energy beams available today or 
will be available in the nearest future: CNGS, NuMi, JPARC (T2K setup) and TeVatron 
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Mat, GeV Mjv, GeV Mjv, GeV 

FIG. 13: Branching ratios of semileptonic decays B — > DlNj (black lines), DglNi (dark gray 

lines) and B^ — > r^clNi (light gray lines), / = e,n,T (solid, long dashed and short dashed lines) 
as functions of heavy neutrino mass Mjv in models: a) I, b) II, c) III. In a phenomenologically 
viable model and heavy neutrino mass within M-j^ < Mjv ^ Mb, the branching ratios are confined 
between corresponding thin and thick lines which show upper and lower limits on from Fig. 4b, 
respectively; form factors for i^^-meson decays are taken to be equal to the form factors for S-meson 
decays from Ref. [71], for Bc-meson decays we adopted form factors from Ref. [72]. 



(NuTeV setup). The relevant parameters of these beams are presented in Table I. To 



Experiment 


E, GeV 


NpoT, 10^9 


Mpp [68] 


Xs [66] 


Xc [67] 


Xb [67] 


(pf), GeV 


(pF), GeV 


(pf ), GeV 


CNGS [73] 


400 


4.5 


13 


1/7 


0.45 • lO"-'^ 


3 • 10-s 


44 


58 


58 


NuMi [74] 


120 


5 


11 


1/7 


1 • 10-^ 


10-10 


24 


24 


24 


T2K [75] 


50 


100 


7 


1/7 


1 • 10-^ 


10-12 


8.5 


10 


10 


NuTeV [76] 


800 


1 


15 


1/7 


1 • 10-3 


2 • lO-'^ 


68 


82 


82 



TABLE I: Adopted values of relevant for heavy neutrino production parameters of several experi- 
ments. 
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FIG. 14: Branching ratios of semileptonic decays B D*lNi (black lines), Bs D*lNi (dark 
gray lines) and Be — J/il^lNj (light gray lines), I = e,fi,T (solid, short dashed and long dashed 
lines) as functions of heavy neutrino mass M/v in models: a) I, b) II, c) III. In a phenomenologically 
viable model and heavy neutrino mass within M^^ ^ Mj^ ^ M^, the branching ratios are confined 
between corresponding thin and thick lines which show upper and lower limits on f/^ from Fig. 4b, 
respectively; form factors for Sg-meson decays are taken to be equal to the form factors for 5-meson 
decays from Ref. [71], for Bc-meson decays we adopted form factors from Ref. [72]. 

estimate the mean longitudinal momenta {ph,l) of D- and S-mesons we make use of the 
parameterization 



with c = 7.7 for E = 800 GeV [77, 78], c = 4.9 for E = 400 GeV [79] and c = 3 for 
E — 120 GeV and = 50 GeV as an factor-of-two estimate. In case of kaons we use the 
estimate 



As we show in Section V, the dominant contribution to the total neutrino production in 
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FIG. 15: Branching ratios of semileptonic decays B — > ttINj (black lines), B — > plNj (dark gray 
lines) and Bg K*lNi (light gray lines), I = e,fi,T (solid, long dashed and short dashed lines) 
as functions of heavy neutrino mass Mjv in models: a) I, b) II, c) Ill.In a phenomenologically 
viable model and heavy neutrino mass within M-j^ < Mjv ^ Md, the branching ratios are confined 
between corresponding thin and thick lines which show upper and lower limits on C/^ from Fig. 4b, 
respectively; form factors are taken from Refs. [70, 71]. 

collisions come mostly from two-body hadron decays. Thus, with neutrino longitudinal mo- 
mentum uniformly distributed in hadron rest frame one gets for average neutrino momentum 
in laboratory frame 



ipN,L)H^:^{pH,L)-[l + ^J . 



If neutrino decay length exceeds detector length AZ, the total number of neutrino decays 
inside the fiducial volume is 



decays 
N 



Al 



M 



N 



tn ^ {pn,l)h 



-N 
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FIG. 16: Inclusive heavy lepton production by strange (black lines), charm (dark gray lines) and 
beauty (light gray lines) hadrons in models: a) I, b) II and c) III; within i^MSM the interesting 

rates arc between corresponding thin and thick lines which show upper and lower limits on 
from Fig. 4b, respectively. 

with being a relative contribution of a given hadron H to total neutrino production, 

h^ Nh{E)-Bt{H^N...) 

where the number of produced hadrons of a type H is estimated as 

Nh{E) = Npot{E) ■ M„{E) ■ xq{E) -BriQ^H) . 
Finally we obtain for the total number of neutrino decays inside the detector 



For the four available beams with parameters presented in Table I and AZ — 5m the 
quantitative predictions are given in Fig. 17 for the three benchmark models with account of 
all experimental and theoretical constraints; both lower and upper bounds scale with mixing 
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Mtv, GeV Mn, GeV 



FIG. 17: Number of sterile neutrino decays within 5m-length fiducial volume for a) CNGS, b) 
NuMI, c) T2K, d) NuTeV beams as a function of sterile neutrino mass M^. Black, dark gray and 
light gray lines refer to benchmark models I, II and III, respectively; in phenomenologically viable 
models the number of decay events are confined by corresponding thin (upper limits) and thick 
(lower limits) lines. 

as C/^. One has to multiply these numbers by the value of the corresponding branching ratio 
(see plots in Section IV), if interested in a particular neutrino decay mode. 

Note in passing that in these estimates we neglected neutrino beam spreading due to 
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nonzero average transverse momentum {pn,t)- With a detector of width /SIt placed at a 
distance I from a target this is justified if 

Otherwise the predictions for neutrino signal have to be corrected by a suppression factor of 
order Similar suppression factor should be accounted for in case of off-axis detector. 

Likewise, as the next approximation to TV^^*^"^* one has to consider the realistic distributions 
of hadron and neutrino momenta instead of their average values. Finally, in case of real 
neutrino experiment, additional suppression can emerge due to focusing system. For a given 
experiment with fixed geometry and detector not tuned to searches for heavy neutrinos, 
these suppression factors can be obtained after dedicated studies; they can be as large as 
one-two orders of magnitude, depending on the neutrino mass. This studies are beyond the 
scope of this work. 



VII. CONCLUSIONS 



A pair of relatively light and almost degenerate Majorana leptons is an essential ingredient 
of the i/MSM. In this model these particles are responsible for observed pattern of neutrino 
masses and mixings and for baryon asymmetry of the Universe. In the present paper we 
analysed the properties of the singlet fermions of the i/MSM, which can be used for their 
experimental search. In particular, we discussed their production in decays of different 
mesons and in pp collisions. We studied the decays of singlet fermions in a wide range of 
their masses and other parameters, consistent with cosmological considerations. 

In the z/MSM the strength of the coupling of singlet fermions to ordinary leptons is 
bounded both from above by cosmology (baryon asymmetry of the Universe) and from below 
by neutrino oscillation experiments. In addition, a lower bound on the strength of interaction 
comes from Big Bang Nucleosynthesis. We analysed the latter constraint in some detail and 
showed that it is stronger than the one coming from neutrino experiments for masses smaller 
than 1 GeV. 

The only particle physics searches for neutral fermions that were able to enter into the 
cosmologically interesting region of masses and couplings of neutral fermions, described in 
this paper, is the CERN PS191 experiment [46, 47]. It was performed some 20 years ago 
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and since then no improvements of the bounds were made. This experiment, together with 
the BBN considerations, indicates that the masses of neutral leptons should be larger than 
the pion mass, though definite exclusion of the region below 7r-meson would require more 
theoretical (BBN analysis) and experimental work. Quite interestingly, with an order of 
magnitude improvement of the bounds on the strength of interaction the whole region of 
masses below kaon mass can be scanned and the neutral leptons could be either ruled out 
or found in this mass range. It looks likely that the significant improvement of the results 
of [46, 47] can be achieved by the reanalysis of the existing data of KLOE and of E949 
experiments, and that the third stage of NA48 would allow to settle down the question. 

The search for the neutral fermions that arc heavier than K-mesons would require ded- 
icated experiments similar to the CERN PS191 experiment, but with higher intensity and 
higher energy proton beams. We argued that the use of the CNGS, NuMI, T2K or NuTeV 
beams can allow to touch the interesting range of mixings for the lepton mass below charm, 
whereas for going above charm much more intensive accelerators would be necessary. 

In conclusion, it is quite possible that the already existing machines can be used for 
the search of the physics beyond the Minimal Standard Model, responsible for neutrino 
oscillations, dark matter and baryon asymmetry of the Universe. Clearly, to study the CP- 
violation in interactions of neutral leptons more statistics would be required, calling for the 
intensity (rather than energy) increase of the proton beams. 
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APPENDIX A: STERILE NEUTRINO DECAYS 



Formulae for most decay rates presented here can be obtained straightforwardly by mak- 
ing use of the formulae for Dirac neutrinos from Ref . [54] ; we checked them and found to be 
correct. Formulae for decay rates into fixed final states are identical in Dirac and Majorana 
cases. Decays N —>■ ti^v and N ^ pv have not been considered there. Our result for decay 
N — >■ 7r°z^ differs from the estimate used in Ref. [55] by an additional phase volume factor 
and by a factor 1/2. Hence, the neutrino life-time used in Ref. [55] to get BBN limits should 
be multiplied by 2/(1 — M^/M^), and corresponding hmit on neutrino mixing angle should 
be divided by this factor. 

Two-body decay modes are 



r (TV ^ n%c) 



r (TV ^ H+1-) = \E^G%\VHfffjM'^ 



16n 



N ■ 



1 - 



Ml 



Ml 



1 + 



\ 



If/a 



327r 



{Mh - Ml] 
Ml 



G%f^,M% ■ 1 



{Mh + MiY 
Ml 



327r 



Ml^ 
Mhj 

Mf 
'Ml 



r(Ar^p+Z-) = ^^G^|KdpM3 



1 - 



Mf 
Ml 



+ 



Ml 
Ml 



1 + 



Mf 



2Ml 



Ml 



\ 



1 - 



{Mp ~ Mif 
Ml 



{M. + MiY 



r (JV ^ p%^) 



Ml 
2 t 

M2 



N 



Ml 


1 t 




Ml 


{- 




Ml 



where Gp is Fermi couphng constant, = I-S/tt, /,,' = —0.45/^, Qp — 0.102 GeV^ [68]; 
hereafter and for CKM matrix elements we use values from Ref. [68] , while for meson decay 
constants we used most recent values from Refs. [68, 80]. 



11 


7T+ 


K+ 


D+ 


D, 


D+ 




Be 


fH, MeV 


130 


159.8 


222.6 


280.1 


190 


230 


480 


Vh 


Vud 




Vcd 


Vcs 


Vub 


Vus 


Vcb 
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Three body decay modes read 



r(A^-W.%)- 192.3 

r (n .jp-^) = ^ 



\Ua\'^{l-8xf + 8xf-xf-12xfl0gxf) , XI 



max 



[Mi„, Mi^ 



M 



N 



1927r3 ''^"l 



with 



and 



(Ci • (1 - S^p) + Cs ■ 6af3) (1 - Uxf - 2x1 - 12a^r) \/l - 4a;; 



+ 12x1 {4 - 1) + 4 (Ca • (1 - 5a/3) + C4 ■ Sai3) (^xf (2 + lOxf - 12xf) ^l-4a;2 
+ 6a;f (1 - 2xf + 2xf) L 



L = log 



(1 + ^ 



1 - 4a;; 



a;; 



Ml 
Mn 



Ci = -{1 - 4sin2 + Ssin* 0^) , 
= 1 (1 + 4sin2 9^ + Ssin^ , 



Co 



■ sin 9„„ (2 sin 



C4 = - sin^ (2 sin^ 6*^ + 1) . 



The Majorana neutrino total decay rate is a sum of all rates presented above multiplied 
by a factor of 2, which accounts for charge-conjugated decay modes. 

APPENDIX B: DECAYS INTO STERILE NEUTRINO 



Differential branching ratio of pseudoscalar meson leptonic decays into sterile neutrinos 

reads 



dBi {H+ l+N) 



Ml 

Mh 



Mlf Ml 
Ml^MlV Ml 



,.M%_Ml 
Mjj Ml 



Mjj - Mf + Ml 



N 



2Mh 



(Bl) 



where th is the meson hfe-time [68]. 
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For differential branching ratios of pseudoscalar meson semileptonic decays one has 

dBr{H^H'l+N) _ \Vhh'?GI f ,2 ( .2 .2. {,2(^2 ,^2^ (^2 ^2^2^ 

+ 2/+((?2)/_ (g2) (2M^ - 2Mfj, - AEnMh - Mf + + q^) + Mf {AEnMh + Mf - Mf, - g^)) 
/2 (^{4EnMk + Mf -Ml- q^) (2M| - 2M^ - AEnMk - Mf + + q^) 

- {2Ml + 2Ml - q') (g2 - - Mf) 



(B2) 

where = {pi + Pn)'^ is momentum of leptonic pair, Vhh' is corresponding entry of CKM 
matrix and f+{q^), f-{q^) are dimensionless hadronic form factors [68] can be found in 
hterature. 

For three-body decays into vector mesons V one obtains 

■ \Ua\ ■ oo,.^A^__ X dq [ — [q -M^-Mi +u 



dEj, - 327r3MH y 2 ^ ' M^ 



"■v 



+ 4{Ml, + Mf) {q' Ml + Mf) qA + 2/|Af^ ( ^ - q') ( M^ + Mf - q' + co' 



+2/3/5 (M^c.^ + (n^ - c^) ) - q'^ + 2hh [q' (2a;2 - O^) + 0^ (M^ - Mf)) 

. f2h f 2^' 



I (M^ - Mf) + ^Mf + 2 {Ml- Mf)' - 2q' (M^ + )) 

+/2/3 (fi'^'^^^ + 2a;^ (Mf - Ml) + {Ml - Mf - q')^ 
+ff {n^ {q^ -Ml + Mf) - 2M^ [q^ - {Ml - Mf)') + 2uj^fl^ {Ml - q' - Mf) + 2u;V)^ , 

(B3) 

where = Mjj - Mfr + Ml - Mf - 2MhEn and Vt^ = Ml - Mf, - q^; form factors 
fi{q^) can be expressed via standard axial form factors AQ{(f') , Ai{q^) , A2{q^) and vector 
form factor V{q^) as 

h = ^ , h = {MH + Mv)-Ar, h = - 



Mff + My' ^ " " Mh + Mv' 

U = (My (2Ao - - A2) + Mff (A2 - ^1)) • ^ , h = h + h 

g2 



which can be found in literature. 
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For two-body decays of r-lepton into heavy neutrino and meson we obtain 

— — ) \ — — J''r" — 



dBv (t ^ pN) _ 



1-^ + 



M2 



M2 



2M^ 



where r,- is r-lepton hfe-time. For three-body decays of r-lepton one has 

rf^;^ -^^-^^^^^-^^V + M2 -^M^J - M2 + M2 - 2i;^Mj V^iv-M^, 



,f2 

-'at 



' (M^-j;jv)' El -Ml 

Mr 2,Mr 



^^Ml-2ENMr 



Note that omitted here charge-conjugated processes also contribute to Majorana neutrino 
production. 
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